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We show how the recently discovered solvable rational extensions of Harmonic Oscillator and 
Morse potentials can be constructed in a direct and systematic way, without the need of supersym- 
metry, shape invariance, Darboux-Crum and Darboux-Backlund transformations. 
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,.1^". I. INTRODUCTION 

T— I \ 

' It is fair to say that in the last three years some of the most interesting developments in mathematical physics have 
" been the discoveries of new types of orthogonal polynomials, called the exceptional orthogonal polynomials, and the 
O ' quantal systems related to them [1-18]. Unlike the classical orthogonal polynomials, these new polynomials have the 
2^ remarkable properties that they still form complete sets with respect to some positive-definite measure, although they 
start with degree £>1 polynomials instead of a constant. 
' Two families of such polynomials, namely, the Laguerre- and Jacobi-type Xi polynomials, corresponding to £ = 1, 
fS| , were first proposed by Gomez-UUate et al. in [l|, within the Sturm-Lioville theory, as solutions of second-order 
eigenvalue equations with rational coefficients. The results in [l[ were reformulated in the framework of quantum 
mechanics in Q , and in supersymmetric quantum mechanics using superpotential in Q . These quantal systems turn 

• out to be rationally extended systems of the traditional ones which are related to the classical orthogonal polynomials. 
I I The most general Xi exceptional polynomials, valid for all integral t = 1,2,..., were discovered by Odake and Sasaki 

1-^, . [J (the case of £ = 2 was also discussed in [3]). Later, in Q equivalent but much simpler looking forms of the 
' Laguerre- and Jacobi-type Xg polynomials were presented. Such forms facilitate an in-depth study of some irnportant 
: properties of the Xg polynomials. Very recently, such systems have been generalized to multi- indexed cases (igI. [l7|. 
I— —I' Other rational extensions of solvable systems, which are not related to the exceptional polynomials, are possible 
\ [T^ - [23| . In these systems, the polynomial part of the wave functions start with degree zero. One of the simplest 
' example of such systems was discussed in [20j, which was later shown to be a certain supersymmetric partner of the 
harmonic oscillator in [2l|. Extending the superpotential scheme of [2l[, Gandati and Berard were able to generate 
an infinite set of solvable rational extension for translationally shape- invariant potentials of the second category j23 | . 
\^ ' More recently, rational extensions of the Morse and Kepler-Coulomb potentials have also been obtained by means of 
Qf^ , the Darboux-Backlund transformation in [23j . 

jy-^ From the viewpoint of the generalized Crum's theorem, solvable rationally extended systems related to the ex- 

, ceptional polynomials are obtainable from the corresponding ordinary systems, which are related to the classical 
^— I ■ orthogonal polynomials, by deleting the lowest energy levels including the ground states [13. The rational extensions 
T— I ' of the harmonic and isotonic oscillators considered in [20l- [23 can be obtained in the same way, with the exception 

• that the ground states were not deleted (see Appendix A of |24l|). 
So far most of the methods employed to generate solvable rational extensions of ordinary systems with or without 

the exceptional polynomials have invoked in one way or another the ideas of shape invariance and/or the related 
Darboux-Crum transformation (supersynimetry) . This requires an exactly solvable ordinary system, such as the 
harmonic oscillator, to be known in the first place, and the superpotential associated with such system is modified 
for the extension. 

In isj we have proposed a simple constructive procedure to generate the exceptional orthogonal polynomials without 
the need of shape invariance and Darboux-Crum transformation. Thus in our work an exactly solvable ordinary system 
and its associated superpotential need not be assumed a priori as in the other works. The superpotential, as well 
as the potential, the eigenfunctions and eigenvalues of the new system are all derived from first principle in our 
method. To distinguish the different roles the superpotential play in our approach and in those employing Darboux- 
Crum transformation or supersymmetry, we prefer to call the superpotential "prepotential" , and our procedure the 
"prepotential approach" . 

It is the purpose of this paper to demonstrate that the solvable rational extensions of the harmonic oscillator given 
in (20r 22] and the Morse potential in [23j can also be generated very simply in the prepotential approach, without 
the need of supersymmetry, shape invariance, Darboux-Crum and Darboux-Backlund transformations. These two 
systems are in the same class as the rationally extended Jacobi system discussed in Sect. 4.4 of [Ts] . 
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II. PREPOTENTIAL APPROACH 



The main ideas of the prepotential approach are summarized here. We refer the reader to [18[ for the details of the 
procedure. We adopt the unit system in which H and the mass m of the particle are such that H = 2m = 1. 

Consider a wave function (j){x) which is written in terms of a function W{x) as (f>{x) = exp{W{x)). The function 
W{x) is assumed to have the form 

W{x,r,) = Wo{x)-lnav)+'^Mv)- (1) 

Here r]{x) is a function of x which we shall choose to be one of the sinusoidal coordinates, i.e., coordinates such that 
r]{x)'^, where the dot denotes derivative with respect to x, is at most quadratic in rj. 

The functions Wo{x), ^(77) and p{ri) are functions to be determined later. We shall assume ^(77) to be a polynomial 
in ry. The wave function is 

pWoix) 

m = ^Piv)- (2) 

Operating on (j>{x) by the operator —dp/dx^ results in a Schrddinger equation = 0, where Ti. = —dP/dx^ + V^ V = 
+ W . For simplicity of presentation, we shall often leave out the independent variable of a function if no confusion 
arises. 

Since W{x) determines the potential V", it is therefore called the prepotential. To make V exactly solvable, we 
demand that: (1) is a regular function of x, (2) the function ^(77) has no zeros in the the ordinary (or physical) 
domain of 77(2;), and (3) the function ^(77) does not appear in V . 

For ^(77) = 1, the prepotential approach can generate exactly and quasi-exactly solvable systems associated with the 
classical orthogonal polynomials [23. The presence of ^ in the denominators of (j){x) and V{x) thus gives a rational 
extension, or deformation, of the traditional system. We therefore call ^ (77) the deforming function. 

If the factor exp{W[){x)) / ^{rf) in Eq. ([2]) is normalizable, then p{rj) = constant (in this case we shall take p{rj) = 1 
for simplicity) is admissible. This gives the ground state 

M^) = -j^- (3) 

However, if exp{Wo{x)) /S,{r]) is non-normalizable, then (I)q{x) cannot be the ground state. In this case, the ground 
state, like all the excited states, must involve non-trivial 7^(77) 7^ 1. Typically it is in such situation that the exceptional 
orthogonal polynomials arise jlSi] . The cases considered in [23| , which we shall rederived by means of the prepotential 
approach in this paper, are such that exp(Wo(a;))/^(77) is normalizable and thus (poix) is the ground state. 
Following the procedure in [l^, we assume ^(77) to satisfy the equation 

c2{v)e+ci{7i)e+£m^o. (4) 

Here the prime denotes derivative with respective to 77. We choose 02(77) — ±77^, and ci is determined by 

'1 d 



(5) 



where Q{ri) = WoV- The function £"(77) was taken to be a real constant in [3], but here we allow the possibility that 
£{r]) may be a function of rj. Nonetheless, it turns out that all formulae presented in [l^ remain intact. 

By matching Eq. Q with the (confluent) hypergeometric equation, one determines f , Q{ri) and ^(77). Integrating 
Q{x) = Wi)f] then gives the prepotential Wo{x): 

dn^y (6) 

The function ^(77) is then given by a linear combination of ^ and 



(7) 
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Here the two functions F{ri) and G{ri) are determined by 

F{7j) = C2{v)V{v), 

G{r^) = (ci-4)V-C2V', 

with the function V{ri) satisfying 

C2V" + (24 - ci) V' + k -c[+£±£ 



V = 0. 



(8) 
(9) 



(10) 



By matching Eq. (jlOp with the (confluent) hypergeometric equation, one determines V, and thus F{ri)^ G{7]), p{r]) 
and £. 

Once all the relevant functions and parameters are determined, we would have constructed an exactly solvable 
quantal system Hcj) = £(j) defined hy H = jdx^ + V{x), with the wave function ([2]) and the potential 



2?)^ ( I ) - (2W^o?) + ^7) ±ci 



±£. 



(11) 



Lastly, we note that the functions (here we add a subscript to distinguish p corresponding to a particular 
eigenvalue £) are orthogonal , i.e.. 



W'^(x{ri)) 



(12) 



in the 77-space with the weight function 



exp 



dx Wo 



(13) 



III. HARMONIC OSCILLATOR 

Let us choose rj{x) = a; G (—00, 00). Then 77^ = 1. For C2 and ci, we take the upper signs in ci and C2 (it turns out 
that the lower signs give the same model) . Thus C2(ry) = 1 and ci — —2Q{ri). 
Eq. (U]) becomes 

e - 2Qiv)e + £^ = 0. (14) 
Comparing Eq. (I14[) with the Hermite equation 

H';ir,)-2r^H',ir,) + 2emir,) = 0, £ = 0,1,2,..., (15) 
where Hiirf) is the Hermite polynomial, we would have 

av)=Uv\o^)^H,[i^). f~ = 2£, Q{V) = V- (16) 

But this choice is not viable, as (77; a) — Hiij]) has zeros in the ordinary domain (—00, 00), which we want to avoid. 
A simple way to solve this is to make the zeros of Hermite polynomials lie on the imaginary axis. This is achieved if 
we set rj ir] in Eq. (fTSjl . giving 

H'/{ir,) + 2rjH^iifj)-2£Heiifi) = 0, £ = 0,1,2,.... (17) 

Matching Eq. (HH with gives 

av)=Uma)^miv), £ = -2i, Q{v) = -v- (18) 

One notes that the Hermite polynomials are odd functions in 77 for odd £. So in this case has a zero at 77 = 0. 
Further study of this case reveals that the wave functions are not normalizable. So here we shall only consider the 
case with even £ — 2m{m =1,2,.. .). By Eq. ([6]), the form of Q{rj) leads to 

W,{x) = (19) 
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We shall ignore the constant of integration as it can be absorbed into the normalization constant. As noted in Sect. II, 
in this case ^(77) = 1 is admissible, as exp{Wo{x))/£^{r]) is normalizable. So the energy and eigenfunction of the ground 
state of this system are So = and 4>o{x) = exp(Wo(a:))/^(77). Below we determine the energies and eigenfunctions of 
the excited states. 

With the solutions in Eq. , Eq. (fTO)) becomes 

(20) 



V" - 2riV' + " 2£ - 2) V = 0. 
Comparing Eqs. (|20|l and (fTSjl (with £ replaced by n = 0, 1, 2, . . .), we get 

V(?7) = i?„(r?), £ = £e^n^2in + e+l), i = 2m. 
From Eqs. ([8]) and ([9|), one eventually obtains 



(21) 



(22) 



Use has been made of the identity = 2-qHn — Hn+i in obtaining the last line in Eq. p2|) . We note that pe,n{'il) is 
a polynomial of degree £ + n + I. 

By Eq. ([T2|) . one finds that Pf,n(?7; a)'s are orthogonal in the sense 



df] — 2-p^,„(ry;a)p£,fe(77;a) cx Snk- 



(23) 



The exactly solvable potential is given by Eq. (fTTj) with 14^0(2^) and ^e{ri; a) given by Eqs. (|19p and ([18]), respectively. 
Explicitly, the potential is 



V{x) =x^-l + 2^ 



2r, 



2£. 



The complete eigenfunctions and energies are 



(24) 



(j)o(x\a) cx 



(j)i^n{x]a) cx 



e 2 



e 2 



it 



0, 



■pLn{ri{x)\a), £Ln = 2{n + £+l), n = 0,l,2. 



Using the identity for £ = 2m [26|, l27| , i.e.. 



H2„^{iii) = (-l)'"22™m!LL ^^(-?7'), 



where L^"-* (77) is the Laguerre polynomial, and the identity 

d 



drj 



we can reduce Eq. (^5)) to 

'/'2m, n 



e 2 



iL^^(-r/2) 



iiL ^^(-r;2)ff„+i(r;) + ry4f!i(-r;2)ff„(r;) 



(25) 
(26) 

(27) 

(28) 

(29) 



This result is identical with that given in [23|. 
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IV. MORSE POTENTIAL 



Now we consider rational extension of the Morse potential. It turns out that in this case £ cannot be a constant. 
Let us choose ri{x) = G (0,cxd), with rj^ = rf. For definiteness we shall take the upper signs for C2 and ci, as 
the lower signs lead to the same results. So we have €2(77) = and Ci = (77 — 2Q{rf)). 
Equation determining ^ is 



In order to link Eq. (j30p with the Laguerre equation 



ryL"^") + (a + 1 - r/) 4^"^ + ^Lf ^ = 0, £ = 0,1,2,.. 



(30) 



(31) 



we rewrite Eq. (|30]) as 



1 - 2 



1] 



(32) 



Directly matching this equation with Eq. ([3T|) will leads to unnormalizable wave functions. So instead we set 77 —rj 
in Eq. This gives 



^r(-^)+ 1 + 2 



,Q(-'7) 



)e'(-.) + ^e(-.) 



-7;) = 0. 



Comparing this equation with Eq. ([5l]) leads to 



or equivalently. 



^e{v;a) = £(77)^-^77, (3(?7) = -^77 - ^J?, 



The form of Q{ri) leads to 



Wo(x) 



a 77 

- — In 77 . 

2 ' 2 



(33) 



(34) 



(35) 



(36) 



According to the Kienast-Lawton-Hahn's Theorem [2a,[23|, the deforming function ^^(77) will have no positive zeros 
in (0,00) if: (i) —2k — 1 < a < —2k with —£< a < —1, or (ii) £ is even with a < —i. Again, in this case, 
^(77) = 1 is admissible. Thus the energy and eigenfunction of the ground state of this system are £q = and 
(poi^) — exp(Wo(x))/^(r/). We now determine the energies and eigenfunctions of the excited states. 
With the solutions in Eq. ([551) . Eq. pUj) becomes 



V" + (-a + 3 - 77) V + 



£-a + l 



-{£ + 2) 



V = 0. 



(37) 



In order that £ be dependent on n, we try V — TfU{'q) where 7 is a real parameter and U{'q) a function of 77. From 
Eq. ^ we get 



77t/" + (27 - a + 3 - 77) [/' + 



f - a + 1 + 7(7 - a + 2) 



(7 + ^ + 2) 



[/ = 0. 



(38) 



Matching this equation with Eq. ([31]) . we have {n = 0, 1, 2, . . 



7 = -(77 + ^+2), 

£i,n = a - 1 - 7(7 - a + 2) = a - 1 - (77 + + 2){n + £ + a), 

P = 27- a + 2 = -a - 2(71 + ^+ 1), 

f/„(77) = L^(77), (3 > -1. 



(39) 
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Putting all these results into F{i]) and G{r]) gives 



(40) 



Pi^rj] a) is a polynomial of degree £ + + 1. It is also easy to check that pe,n{ill Q^)'s are orthogonal with respect to 
the weight function 



The exactly solvable potential is given by 



V{x) = + ha ~U- l)e-" + ^ + 2^' 

A 2 4 



-2x 



ae 



The complete eigenfunctions are 



00 (a;, a) oc 



e 2 7y 2 



e 2 7^ 2 

(f>e,n{x;a) oc p£,„(?7;a), 



(41) 

(42) 

(43) 
(44) 



where pf^„(?7;a)'s are given in (|40l) . The corresponding eigen-energies are £0 = and in (|39|) . respectively. For 
the wave functions to be regular at x = 0, one must have n < —a/2 — £—1. This means the system admits only finite 
number of bound states. 

With the help of the identities and 



riLtX'\ri)-aLf\{i^) = -lL^r'\^)^ 

one can recast Pi,n{'ri\ Q^) in l|40|) into 

Pi Am «) = - [(" + ^) Lf\{-r,)L\^P\^) + (n + l)Lf V^)4'i W ■ 
This expression is exactly the same as that given in (2^ in the case i — 2m, with the identification 

a = -2(a + ^+l), 
n ^ k, 

13 = -a-2(n + ^+l) = 2(a-fc). 



(45) 
(46) 

(47) 



(48) 
(49) 



V. SUMMARY 



We have shown how the recently discovered solvable rational extensions of Harmonic Oscillator and Morse po- 
tentials can be constructed in a direct and systematic way, without the need of supersymmetry, shape invariance, 
Darboux-Crum and Darboux-Backlund transformations. In our approach, the prepotential, the deforming function, 
the potential, the eigenfunctions and ei gen values are all derived within the same framework. 

With the results given here and in [l8|, rational extensions of all well-known one-dimensional solvable quantal 
systems based on sinusoidal coordinates have been generated by the prepotential approach. One would like to apply 
the same approach to find rational extensions of the other solvable models based on non-sinusoidal coordinates, 
following the work of the third paper in [25| . Unfortunately, such way of rational extensions only lead to quasi-exactly 
solvable systems, because the energy quantum number n will appear in the x-dependent terms in V{x). A non-trivial 
generalization of the present approach may be in order, which we hope to report in the near future. 
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